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Abstract 

Atomic and electronic structures of TiO2(110) surface with possible adsorp- 
tional, substitutional and interstitial Au or Pt elemental impurities at full and 
one-sixth monolayer concentrations were investigated by density functional the- 
ory calculations using the projector augmented wave approach within the plane 
wave method. Relative thermodynamic stabilities of such phases have been dis- 
cussed by means of their surface free energies. Our results suggest that tunable 
photocatalytic activity can be achieved on Pt atom admixed rutile (110) surface 
at low coverages. 
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1. Introduction 

Titanium dioxide (Ti02) attracts an increased interest as a catalyst sup- 
port. Gold nanoparticles on titania were shown to exhibit remarkable catalytic 
activity [l[. Another important mctal-on-oxide system is platinum incorpo- 
rated Ti02 which finds applications as catalysts and gas sensors 01. Pt/Ti02 
system has drawn further attention because of its photocatalytic activity to- 
wards water decomposition In recent studies, these noble metals supported 
on Ti02 has shown to perform highly efficient catalysis under solar light irradi- 
ation 0,1,0,0]. 

As a wide-gap semiconductor, the rutile TiO2(110) surface is considered to 
be the generic model system for oxide surfaces. The (110) surface of rutile struc- 
ture has the lowest surface energy among the other facets [ii,[i3,[ii|,[i2. It IS 



reducible by surface oxygen vacancy creation or metal incorporation which at- 
tracts a great deal of interest for fundamental study of photo- and heterogeneous 



mi 



tracts a great aeai oi mterest tor lunaamcntal study or pnoto- and netero 
catalysis [13, Ei, 15, 16|, functional ultrathin films [l3,[l3 ^^"^ dielectrics 
Understanding of the properties of metal-metal oxide interface can provide im- 
portant insights into the applications of real catalysts. 
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Bonding of gold on titania has been studied to shed hght on the effect of 
Au-su ppo rt adhesion on its catalytic behavior both exp er iment ally 12 1[ 1 2 2| . 1 2 3l 
m m M [23 and theoretically [H, [H, M, M, M, S, M, M, S, M, M, S Ej- 



Similarly, the adsorption properties of Pt haye also been investigated b y many 
experiments 0, IH 111 [H, 13 li^, lii] and by a few theoretical studies [47l [32. 
ksj . Although the nobel metal-enhanced catalytic activity of rutile Ti02 has 
been reported by experiments, the electronic properties of such systems as the 
originating factor are still not well known. 

Moreover, the composition of nobel metals like Au and Pt to Ti02 support 
can be interstitial or substitutional as well as being adsorptional. Indeed, ex- 
periments showed that Pt atoms can thermally diffuse into Ti02 lattice under 
oxidizing atmosphere [4§|. Futhermore, these diffused Pt atoms can substitute 
Ti"*"*" when oxidized to Pt^^ or they form interstitials inside. Such metal impu- 
rities are known to greatly influence the electronic and catalytic properties of 
the combined system. 

In this paper, we studied the structural and electronic properties of adsorp- 
tional, substitutional, and interstitial Au or Pt impurities on rutile TiO2(110) 
surface at high (1 ML) and at low (1/6 ML) concentrations. Relative ther- 
modynamic stabilities of such impurity phases have also been discussed. Our 
primary aim is to elucidate the effect of Au(Pt) incorporation on the electronic 
structure of the titania support at the fundamental level. In this sense, we give 
emphasis on the electronic behavior - especially on the photocatalytic activity 
- as a result of the bonding characteristics of such an incorporation rather than 
investigating the formation of such impurities. 



2. Method 



The calculations have been performed by the density functional theory (DFT) 
implementation of the VASP [50[ code. Exchange-correlation energy has been 
approximated by the gradient corrected Pcrdcw-Burke-Ernzerhof (PB E96) func- 
tional 5l|. We used projector augmented waves (PAW) approach [H, ^ with 



a plane-wave basis up to a cutoff of 400 cV. 

We considered the bulk terminated (1x1) and (3x2) supercells as the slab 
models for high (1 ML) and low (1/6 ML) metal coverages, respectively. Bulk 
terminated rutile TiO2(110)-(lxl) and TiO2(110)-(3x2) surfaces have been 
modeled by a symmetric slab of 7 Ti02 trilayers separated by ^^15 A of vacuum 
region. Each trilayer consists of a central 0-Ti-O plane and 2 oxygen atoms 
placed symmetrically above and below this plane. 

For the geometry optimization calculations, the Brillouin Zones of (1x1) 
and (3x2) supercells were sampled with a 8x5x1 and 2x2x1 fc-point meshes, 
respectively. In all calculations, the full relaxation has been performed using 
conjugate-gradient algorithm based on the reduction of the HcUmann-Feynman 
forces on each atom to less than 0.01 cVA~^. We used much denser grids for 
the computations of band structures and densities of states (DOS). 
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We calculated the binding energy of adsorbate, M, by 



Eb = {Eu/TiO^ ^ £'Ti02 ^ 2i?M)/2 

where -Eivi/TiOj is the total energy of the computation cell involving the slab and 
the atomic impurities, M, i?Ti02 is that of the defect-free stoichiometric slab, 
and Em is the energy of an isolated M (Au, Pt) atom calculated in its electronic 
ground state. 

The formation energies of defects in the form of metal atom impurities 
on (1 X 1) and (3x2) surfaces has been calculated (as previously described in 



detail[5J,l55|) by 

Ef = (-BM/Tio^ " -BxiOs - ^m^iM + An/iTi)/2A 

where Am and An are the differences in the number of adsorbate M atoms and 
surface Ti atoms from the reference stoichiometric Ti02 slab, respectively. The 
chemical potentials, /iM, were taken from their reference bulk values of —3.270 
eV for Au and of —6.017 eV for Pt, representing their most stable solid phases 
accessible. In other words, we assume that the surface is in thermodynamic 
equilibrium with cep bulk Au or Pt. By the same token, the surface layer must 
be in equilibrium with the rutile Ti02 bulk which comes into contact with. 
This physical requirement implies /iTi + 2/.to = A*Ti02 restricting, /^xi within 
an interval of allowed values. Chemical potential of Ti can be as high as that 
of its bulk which defines an upper boundary referring to Ti-rich conditions. 
On the other hand, molecular oxygen defines the most stable phase for /^o so 
that = \Eq)^ referring to 0-rich conditions. This choice, therefore, defines 
the minimum value of /xxi through the thermodynamic equilibrium condition, 

MTi + 2^o = MTi02- 



3. Results and Discussion 

Rutile Ti02 has a direct band gap of 3.03 eV [H^ at F corresponding to UV 
optical response. Our calculated value is 1.85 eV because of the well known 
underestimation of GGA functionals due to insufficient cancellation of the self- 
interaction energy. The band gap can greatly be corrected by employing many- 
body perturbative corrections up to first order in the screened Coulomb poten- 
tial, W, called as the GW approximation. A better description of both electronic 
and optical spectra can be obtained by applying Bethe-Salpeter equation (BSE) 
including excitonic effects. Such a quasiparticle treatment for the correlation 
energy starting from the DFT spectrum (Kohn-Sham eigenvalues and wavefunc- 
tions) of Ti02 bulk material has been reported to give scissors-like correction 



to unoccupied states without noticeable change in the band dispersions [57|, |58| . 
Therefore, descriptions of electronic properties based on pure DFT results can 
be made as far as the band structures are concerned. After all, many-body 
perturbative approach for a supercell consisting of TiO2(110)-(3x2) slab with 
7 trilayers and a vacuum region of ~15 A is computationally very expensive, if 
not impossible even at the GW level. 
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Another important factor is that the surface energy tend to converge with 
increasing number of trila yers with an odd-even oscihation similar to those re- 
ported previously [H^ 6^, 61, 62, 63, 64, 65|. After testing the convergence of 



surface energetics as a function of the slab thickness we have chosen the seven 
trilayer model as a compromise between the accuracy and the computational 
cost. The present work not only uses a larger cell (e.g. 252 atoms for (3x2) 
cell) compared to the previous studies but also puts forward the electronic band 
structures of TiO2(110) surface with Au(Pt) atomic impurities at 1 and 1/6 
ML coverages. This allows us to discuss the effects of such a metal incorpora- 
tion on the electronic structures specifically in the gap region to get a better 
understanding. 

3.1. Au(Pt)/TiO2(n0)-(lxl) 

At 1 ML coverage, Au binds to fivefold coordinated surface Ti (Ti5c) with 
a bond length of 3.10 A [dAu-Ti in Tabled] tilted by 14° toward the nearest- 
neighbor (nn) bridging oxygen (02c) as shown in Fig[T] Well-ordered Au mono- 
layers where gold being atop Ti5c on (1x1) surface has also been verified by 
several experiments [Sfl [13, . Single Au adsorbates are separated from each 
other by 2.96 A along [001] direction. Our calculated gas phase Au dimer length 
IS 2.53 A in agreement with Hakkinen et aL's theoretical result of 2.54 A 1661 1. 
Hence, Au adatoms on (1x1) surface form linear metallic chains parallel to the 
oxygen rows where the interatomic distances are slightly longer than the gold 
dimer length. This indicates a strong Au-Au coordination. Projected DOS 
analysis shows that the peak corresponding to the half filled gap state is mainly 
due to Au 5d electrons. It exhibits high dispersion along X-M and X'-F-M 
reflecting the strength of the bonding along ID gold chain. It couples to the 
conduction band (CB) of Ti02 and is almost flat between F and X (reciprocal 
to [110] direction) indicating a weak adsorption. Yang et al. found the binding 
energy of 1 ML Au on the relaxed TiO2(110) to be 1.49 eV at Ti5c site with a 
Au-Ti bond length of 2.66 A through full potential linearized augmented plane- 
wave (FLAPW) calculations [2^. We calculated the binding energy (BE) of an 
Au atom on the stoichiometric (1x1) rutile surface to be 1.38 eV. In fact, this 
BE does not reflect bare metal-substrate interaction, it substantially involves 
interaction of gold with its periodic images on (1x1) cell causing a signiflcant 
increase in adsorption energy. Indeed, when we calculated the BE on (3x2) 
surface it reduces to 0.40 eV where Au adatoms are separated from each other 
by 8.9 A along [001] and by 13.1 A along [110]. Therefore, Au weakly binds 
to the surface. Moreover, calculating the BE's at all possible adsorption sites 
we obtained an almost flat potential energy surface similar to that of Iddir et 
al. [33 . This suggests that gold can diffuse in all directions over the surface in 
agreement with experimental findings [23, 24, 26| . 



Gold substitution for Ti5c at 1 ML causes considerable distortion on the 
surface morphology and on the electronic structure as shown in the second 
row of Fig. [1] The distance between Ti5c and the oxygen beneath it extends 
from 1.83 A to 2.64 A by Au substitution. Metal dopant interacts with four 
nn threefold coordinated basal oxygens (03c) much weaker than Ti does. An 
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r X MX' r M 




Figure 1: Relaxed geometries (top three trilayers) of the Au(Pt)/Ti02 (llO)-(l X 1) systems 
are presented on the left. Their band-gap states along with the corresponding bulk projected 
band structures (shaded areas) and DOS are shown on the right, next to them. Gray and dark 
gray (red in color) balls denote Ti and O atoms, respectively. Some key distances between 
Au(Pt) and (sub)surface atoms are depicted (in A). 
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Tabic 1: Calculated values for the M/TiO2(110) systems (M=Au, Pt); work function and 
Fermi energy relative to bulk valence band top (in eV), as well as M-O and M-Ti distances 



Slab System 




'-^ivi — yj "^ivi — _L 1 


(1x1) clean 


7.22 0.00 




A 

subst 


6.16 0.76 
7.39 -0.40 


2.74 3.10 
1.97 — 


Pt 

subst 


5.00 0.49 
6.78 0.33 


2.02 2.64 
1.96 — 


(3x2) clean 


7.26 0.00 




on 

Au subst 
in 


4.97 1.75 
7.56 -0.05 
5.69 1.80 


2.05 3.59 
2.03 — 
2.18 2.75 


on 

Pt subst 
in 


6.23 0.63 
7.39 0.03 
6.97 0.46 


1.96 2.41 
2.01 — 
1.99 2.70 



explanation might be that the valence of Au, 6s^, compared to the valence of Ti, 
3^*^45^, imposes an electron deficiency. As a result, a half filled Au-03c driven 
state appear within the VB setting the Fermi energy below the bulk projected 
VBM of Ti02. Additionally, an unoccupied impurity state falls in the gap with 
a coupling to the CB at about X. 

Low energy Pt adsorption site is above Ti5c tilted by 24.7° toward the nn 
02c similar to single Au adsorption case but closer to the surface. This site is 
also referred as the hollow site since it is atop the middle point between two 
basal oxygens [s^ . This prediction slightly disagrees with experimental adsorp- 
tion site atop Ti5c, probably due to a difference in the theoretically predicted 
and experimentally measured amounts of charge transfer from Pt to Ti02 (42j . 
Pt adatom causes noticeable lattice distortions up to the second trilayer (Fig.[T]). 
Unlike the gold case, Pt adatom interacts strongly with rutile (110) surface giv- 
ing a BE of 2.79 eV on (1 x 1) cell. It decreases to 2.17 eV on (3x2) cell. Clearly, 
this lowering is much smaller than that of the Au adatom case, since the con- 
tribution from the interaction of Pt with its periodic images is smaller. This is 
because of the fact that Pt-Pt separation along [001] on (1x1) surface is con- 
siderably larger than the Pt dimer length of 2.33 A due to the lattice parameter 
(c=2.96 A) of the bulk Ti02. Moreover, Pt ML reduces the work function of the 
clean surface by 2.22 eV (Table [I]), more than Au does due to larger amount of 
charge transfer from metal to substrate, indicating a stronger binding. Hence, a 
uniform Pt deposition is relatively more probable than an Au adlayer formation. 
In fact, Steinriick et al. reported that Pt deposition happens to cause a uniform 
coverage on the surface at low temperatures (< 160 K). Higher substrate tem- 
peratures lead to Pt islands by increasing Pt diffusion probability [4]J. Pt-02c 
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and Pt-Ti5c bonds bring two impurity states that disperse strongly across the 
band gap. Their flat-hke dispersion along F-X reveals the weakness of Pt~Pt in- 
teraction along [110]. Fermi energy crosses these states leading to metallization. 
Ft induced distortions, localized to surface layer, bring states at and around the 
VBM. An upper lying empty impurity state couples to the CB near and between 
M and X' modifying its edge. 





Figure 2: Top three trilayers of fully relaxed Au/TiO2(110)-(3x2) systems along [110] and 
[001] directions for adsorptional (on), substitutional (subst) and interstitial (in) cases are 
shown. Relevant energy band diagrams are presented with DOS (in the rightmost panel) for 
each of these atomic structures. 



Substitutional Ft on (1 x 1) surface docs not distort the surface layer as much 
as an Au dopant does (at the bottom row of Fig. [1]). The distance between Ft 
and subsurface O is 1.97 A that is slightly extended with respect to Ti-0 bond 
length of 1.83 A. Having fourfold coordination with basal oxygens, Ft dopant 
modifies the VBJVI edge by inducing a gap state that elevates the Fermi energy by 
0.33 eV at F relative to the VBIVI of the clean surface. Another unoccupied state 
appears a '-^0.1 eV above the Fermi energy leading to a narrow-gap (indirect 
between F and X') semiconducting system. 



3.2. Au(Pt)/TiO2(110)-(3x2) 

IVIinimum energy adsorption site for single Au on the (3x2) cell with 21 
layers has been found to be at above bridgi ng oxygen , 02c as shown in Fig. [2] in 
agreement with previous theoretical results j35ll29ll26| . Another preferential site 
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is the the hohow site [3l|, [S^, l37 1 . These two configurations differ insignificantly 
in their total energies and are both experimentally verified [l^l- Campbell et 
al. reported Au monomer adsorption energy to be 0.43 eV by calorimetric 
Vijay et al. found that Au binds to an 02c or to a Ti5c 



measurements 



67| 



atom (with a tilting toward an 02c) weakly by 0.6 eV. We calculated the binding 
energy of Au on the (3x2) cell to be 0.40 eV which indicates a flat-like potential 
energy surface. Indeed, the diffusion barriers were found to be so low that Au 
atoms already diffuse at room temperature (RT) both in the [001] and in the 
[110] directions [H El [H IH . 




Figure 3: Atomic structures of Pt/TiO2(110)-(3x2) systems along [110] and [001] directions 
for adsorptional (on), substitutional (subst) and interstitial (in) cases. Relevant energy band 
diagrams and DOS are presented on the right for each of them. 



On stoichiometric TiO2(110) surface Au is reported to exhibit a quasi-two- 
dimensional growth at low concentrations at low temperatures. [21, 23 1 As the 
deposition rate or temperature increases Au starts to form three dimensional 
islands. Formation of Au clusters on the surface indicates a very weak metal- 
substrate interaction. On the other hand, strong binding can be achieved at 
02c vacancies. For example, Benz et al. deposited Au atoms on titania surface 
forming monatomic Au centers [isl ]. Following experiments by Tong et al. re- 
vealed that Au binding can be broken by a hydroxyl group forcing Au out of 
the adsorption site in the presence of water. 

Electronically, single Au adsorbate causes metallization due to an unpaired 
5d electron. IVIoreover, it gives a state just around the CBIVI and modifies the 
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CB edge. Okazawa et al. determined the work function of TiO2(110) surface 
with low gold coverage to be ~5.3 eV where our value is 4.97 eV. On the other 
hand, their clean surface value of 5.4 eV is nearly 2.2 eV smaller than ours. 
The discrepancy can be addressed to possible existence of oxygen vacancies on 
the sample because such a reduction was shown to cause 2 cV drop in the work 
function by Vogtenhubcr et al. [t^]. Our clean surface value rather agrees with 
theirs of 7.16 cV and also with an experimental result of 6.83 cV as reported in 
Ref. 



III. 



Thermodynamic stability of Au(Pt)/TiO2(110) 
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Figure 4: Calculated formation energies of Au(Pt)/Ti02 (110) systems as a function of the 
chemical potential of Ti. 



When Au substituted for an in-plane Ti5c, it distorts nn oxygens similar to 
(1x1) case, only locally. Furthermore, the overall band structure is very similar 
too. Due to charge deficiency, Fermi energy falls in the VB of Ti02 giving rise to 
metallization. Clearly, in substitutional cases gold also raises the work function 
by 0.17 on (1x1) and by 0.30 eV on (3x2) cells due to valence electron deficit. 

We found that precious metal atom implantation into TiO2(110) subsurface 
interstitial sites at 1 ML concentration is unstable. However, our calculations 
suggest that metal atoms are possible as interstitials in (3x2) cell, probably at 
elevated temperatures and/or under oxidizing conditions in agreement with the 
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experimental findings I43,[Z2|- Relevant surface energy diagrams are presented 



in Fig.|4]where the formation of An in (3x2) phase is about 0.75 eV higher than 
that of the clean surface. Single Au implanted cell appears to be metallic just 
like the previous cases. A flat going impurity state lies '^0.3 cV above the VBM 
of Ti02. 

Adsorbate-induccd modification of the substrate morphology is more pro- 
nounced for Pt on rutile (110) relative to the Au case. For instance, Pt pulls nn 
Ti5c and 02c toward itself at the expense of extending their bonds with the lat- 
tice. This indicates relatively stronger binding with 2.17 eV adsorption energy 
at the hollow site (see Fig. [3]) in agreement with the BE of 2.14 eV reported by 
Iddir et al. for Pt on (3x2) with 4 trilayers. This site differs with a tilting angle 
of ^27° from experimentally assigned adsorption atop Ti5c atom The dis- 
agreement has been addressed to the difference in theoretically predicted and 
experimentally estimated amounts of charge on Pt. Experiments confirm the en- 
hancement of photocatalytic activity of titania surfaces by platinization (45. 44 1. 
This can be partly explained by the availability of defect driven gap states that 
increase transition probabilities. For example, the excess charge localized on 
the Pt adsorbate causes three flat going states in the band gap, which can also 
be seen from their strong DOS peaks. These states lie at 0.04, 0.35, and 0.63 
eV above the VBM at F. Moreover, surface states appear near band edges due 
to local distortions. A significant band gap narrowing of 0.84 eV is predicted 
resulting from upper lying flat-like occupied state (0.63 eV above VBM) and an 
unoccupied surface state falling into the gap from the CB (0.21 eV below CBM 
at F). 

Relaxation of the clean surface cell causes Ti5c atoms to sink. But, sub- 
stituted Pt stays levelled with in-plane oxygens (Fig. [3]). It slightly pulls the 
subsurface 03c up. Pt mediates less distortion to the lattice in comparison 
to the gold case. Moreover, Pt substitution on (3x2) cell is semiconducting 
whereas a single gold dopant leads to metallization. In addition, we found that 
this phase is thermodynamically more stable than the clean surface under oxi- 
dizing conditions (Fig. HJd) . Choi et al. has recently shown that Pt at low doping 
level signiflcantly enhance photocatalytic activity of rutile Ti02 0] ■ This must 
be related to the impurity states. Pt substitute increases the work function by 
~0.13 eV indicating a charge transfer to the lattice that leaves behind empty 
valence levels on Pt atom. Therefore, Pt substitution on (3x2) cell mainly 
brings unoccupied gap states. Two of which are flat and lie just 0.54 and 0.62 
eV above the VBM. A number of impurity induced states fall into gap from the 
CB that show bulk-like dispersion. Together with these states significant band 
gap narrowing mimics the observed increase in photoreactivity. 

Zhang et al. reported that Pt can substitute Ti^"*" under oxidizing condi- 
tions and can also thermally diffuse into Ti02 substrate [i^. Therefore, we 
considered Pt at the interstitial site as shown in the last row of Fig. [31 We cal- 
culated the formation of such an interstitial phase is ~0.65 eV relative to that 
of the bulk terminated surface (Fig. |4}d) in line with the experiment. Although 
it distorts the lattice more than the previous cases, its effect on the energy 
band structure appears to be comparable to the adsorption phase due to similar 
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bonding characteristics. Pt interstitial brings three occupied states with sharp 
DOS peaks implying charge locahzation around Pt. One of these couples to 
the VB of the clean surface. Relative to the VBM, the others lie ~0.12 eV and 
^0.37 eV higher in energy at F. They make their minimum at F and disperse 
up to the flat section along M-X'. Fermi energy is set at the top of the upper 
lying state 0.46 eV above the VBM. As in the case of Pt adsorption, bulk-like 
surface states fall in the gap from the CB by 0.1 eV at F. Resulting band gap 
narrowing of 0.56 eV corresponds to visible optical response. Such point defects 
at the interstitial sites might also be important to get a better understanding 
of the complex behavior of titania. For instance, in order to offer a possible 
explanation for the appearance of Ti3c? defect state observed on (110) surface 
reduced by the loss of bridging oxygen(s), Wendt et al. considered Ti atom at 
interstitial sites that yields a gap state as observed (72j . 

4. Conclusions 

Au atom binds to the stoichiometric rutile (110) surface much weaker than 
Pt does. Theory suggests a three dimensional clustering upon gold deposition 
on the defect-free surface. In this sense, a full Au adlayer is difficult to realize 
whereas Pt coverage is more probable. Au and Pt diffusion into the (1x1) unit 
cell is thcrmodynamically unstable. However, their substitution for Ti5c become 
even more stable than the clean surface at low concentrations under oxygen 
rich conditions. Noble metal incorporated phases on (3x2) cell are found to 
be within the reach of thermal treatment. Therefore, Pt and Au atoms can be 
adsorbed, or doped as substitutes for the fivefold coordinated Ti atoms, and 
implanted into interstitial sites in the lattice, at low concentrations. Formation 
of adsorptional impurities are energetically more favorable than the other two 
phases. 

Both noble metals are expected to promote the catalytic behavior of TiO2(110) 
surface by increasing the reaction probabilities through availability of band-gap 
states. Stoichiometric Ti02 gains metallic character upon single Au atom pres- 
ence due to an unpaired 5c? electron. Moreover, narrowing of the gap towards 
visible region results from impurity driven defect states for Pt which, hence, can 
be used in photocatalysis. The wavelength tuning of photo response of titania 
might be achieved by different types of Pt incorporation at low coverages. 
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